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We compute the diffusion coefficient and the Lyapunov 
exponent for a diffusive intermittent map by means of cycle 
expansion of dynamical zeta functions. The asymptotic power 
law decay of the coefficients of the relevant power series are 
known analytically. This information is used to resum these 
power series into generalized power series around the alge- 
braic branch point whose immediate vicinity determines the 
desired quantities. In particular we consider a realistic situa- 
tion where all orbits with instability up to a certain cutoff are 
known. This implies that only a few of the power series coef- 
ficients are known exactly and a lot of them are only approx- 
imately given. We develop methods to extract information 
from these stability ordered cycle expansions and compute ac- 
curate values for the diffusion coefficient and the Lyapunov 
exponent. The method works successfully all the way up to 
a phase transition of the map, beyond which the diffusion 
coefficient and Lyapunov exponent are both zero. 

I. INTRODUCTION 

Given the task of computing an average, such as a 
Lyapunov exponent or diffusion coefficient of a chaotic 
system, one can take two different approaches. Firstly, 
simulation is usually simple and it provides an answer 
without bothering to understand the topology of the flow, 
but it may suffer from severe convergence problems. 

Secondly, these averages can be extracted from dynam- 
ical zeta functions and their expansions, known as cycle 
expansions. The basic advantage of expanding the av- 
erage over cycles is that the asymptotic limit < — )■ oo is 
already taken from the beginning. Longer cycles provide 
corrections to the results obtained from shorter ones [Q . 

Real success applying zeta functions has so far only 
been demonstrated for quite a restricted class of dynam- 
ical systems [Q-^- The topology of the flow should be 
Markovian — symbolic dynamics may be introduced and 
this symbolic dynamics is of finite subshift type (mean- 
ing that there is only a finite number of forbidden sub- 
strings). In addition the system need to be hyperbolic 
— the stability of cycles is exponentially bounded with 
length. The class of systems complying with these two 



properties is called Axiom-A. This class is far too re- 
stricted to have any major relevance in applications. 

Success in expanding a zeta function depends on its 
analytic structure. Convergence is hampered by singu- 
larities close to the zero being studied. However, if the 
nature of a disturbing singularity is known, one can uti- 
lize this knowledge in a resummation scheme. If the sin- 
gularity is solely due to intermittency the convergence 
problem is thus tamed to a large extent Q . 

To appreciate the relevance of stability ordering of cy- 
cle expansions we imagine a fairly generic system, given 
by some set of differential equations. The problem of 
finding periodic orbits in a systematic way is largely fa- 
cilitated if one has some symbolic dynamics. For a few 
potentials this is possible, for example the x^y^ model [p|, 
the Helium atom Q , the diamagnetic Kepler problem ||^ 
and the anisotropic Kepler problem |^,^. 

For generic flows it is often not clear what Poincare sec- 
tion should be used, and how it should be partitioned to 
generate a symbolic dynamics. Cycles can be detected 
numerically by searching a long trajectory for near re- 
currences. The long trajectory method for finding cycles 
discussed in ||l^,0 preferentially finds the least unstable 
cycles, regardless of their topological length. If you can 
find all cycles with stability Ap less than a certain cutoff 
you can use stability ordered cycle expansions. Stability 
ordering was introduced in |5|Jl^. It has later been stud- 
ied more systematically in |l^4[| . It is much easier to 
implement for a generic dynamical system than the cur- 
vature expansions which rely on finite subshift approxi- 
mations to a given flow. 

A general stability ordered cycle expansion looks like 
X^Sj"" exp(— sZi), where is a monotonically decreas- 
ing sequence but h is not monotonic. In this paper we 
will restrict our attention to maps. (It would then be rele- 
vant to speak of stability truncation rather than stability 
ordering.) The expansion looks like X^t!^'"' '^i-^' where a 
few of the coefficients may be exact whereas the rest are 
only approximate. In particular if the system is inter- 
mittent the number of approximate coefficients greatly 
exceeds the number of exact ones and the main task of 
this paper is to extract the information they carry. More- 
over, we will make use of our a priori knowledge of the 
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power law decay of the exact coefBcients and employ the 
resummation technique of ref |^ to improve convergence. 

We believe that the idea of stability ordering has its 
biggest potential for systems which cannot be described 
by a symbolic dynamics of finite subshift type. But in 
order to identify the problems due only to intermittency, 
we will study a map with complete symbolic dynamics. 



II. THEORY 
A. Averages and zeta functions 

A nice introduction to chaotic averages is found with 
proper references in ; we will take a slightly different 
approach. The reason for this is that the key step in 
[ p5 | assumes that the leading zero of a zeta function is 
isolated. We will try to avoid this assumption by starting 
from an expression for the invariant density in terms of 
periodic orbits. The price we pay is that this formula 
is not rigorously proven for the intermittent systems we 
will study. 

The aim is to compute averages like 



with the zeta function 



{w{x,n)) — / p{x)w{x,n)dz 



(1) 



where p{x) is the invariant density of the ergodic map 
X ^ /(a;)- 

This density can be expressed in terms of periodic or- 
bits via 



p r—l ' ^ ' xi 



(2) 



where r is the number of repetitions of primitive orbit p, 
having period rip, and stabihty Ap = ^^^;^\x=xi with Xi 
being any point along p. 

The weight w{xo,n) is associated with the trajec- 
tory starting at xq and evolving during n iterations 
in such a way that it is multiplicative along the flow: 
w{xo,ni + n2) — w{xo,ni) ■ w(/"^ (a;o), ^12)- As we are 
dealing with maps, it is simply w{xo,n) — wIxq,!) ■ 
w{f{xo), 1) ■w{p{xo), 1) . . . w{f"-^^, 1). The phase space 
average of •w{xo,n) may now be expanded in terms of 
periodic orbits as 



lim (w;(a;o,n)} = lim E^pE^pTXTT ' 

p r-1 I VI 



(3) 



Wp is the weight along with cycle p. Zeta functions are in- 
troduced by observing that the average (|^) may be writ- 
ten as 

lim (w(xo,n)) = lim — / z^'"^ —log Q'^ {z)dz , 
rwoo zm J (J az 

(4) 



i/c4.)=n(i--.^ 



(5) 



C is a small contour encircling the origin in clockwise 
direction. Eq. may be verified by inserting the zeta 
function (^) and let the integral pick up the residues from 
z = 0. The result can be recast into a sum over residues 
outside C, that is, it may be related to the analytic struc- 
ture of the zeta function. 

The Lyapunov exponent can be expressed in terms of 
a generating function 

A= lim -(log|A(a;o,n)|) = hm --^(A(xo,n)'^) l^^o 

n^oo n n^oc n dp 

One therefore introduce the multiplicative weight 

w{xo,n) = A{xo,n)'^ , (6) 

One can now express the Lyapunov exponent in terms 
of the associated zeta function 

For a diffusive map / : R R, the diffusion coefficient 
can also expressed in terms of a generating function 

1 



lim _((/"(io)-xo)') 

n^oo Zn 

rwoo Zn ap'^ 



(8) 



motivating the introduction of the weight 

wi5(a;o,i) ^e''^^"^*""-*") . (9) 

If f{x + nL) — f{x) + nL where i G / (/ is some 
interval of length L) then the map can be reduced to a 
map / : / I— > / on the elementary cell. 

This may be expressed in terms of a zeta function with 
the weight Wp along cycle p on the elementary cell given 
by 



where 



(10) 



(11) 



is the corresponding drift in the full system. 

The diffusion coefficient may now be expressed in terms 
of the associated zeta function 



(12) 
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In both and ( |T^ ) the asymptotic behavior will deter- 
mined by the leading singularity of the integrand. Since 
the integrand is evaluated at /3 = 0, the singularity is lo- 
cated at z = 1. If this singularity is isolated the asymp- 
totic result is obtained by simply integrating around it. 
For the intermittent system we are going to consider 
there is a complication. The singularity is not isolated. 
The zeta function has a branch cut along Re(z) > 1 and 
lm{z) = 0. To extract the asymptotic behavior of these 
integrals we need to integrate around this cut. 

Let us return to the Lyapunov exponent and assume 
that 



B. Expanding zeta functions 

For the Lyapunov exponent calculation we expand the 
zeta function 

= n (i - ^) 

/ ^ z"plog|A| , 



l/Cx{z,l3)^[a,{l-z) + 0{{l-zr)] 

+/3[&o + 0(1 - z)] + 0(/32) (13) 

where 7 > 1. This particular assumption will be moti- 
vated later in this paper. We need to evaluate 
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(1 



dp ds 



ai 



(14) 



where we have changed variable to s = 1 — z. Fq is 
a contour encircling the negative real s-axis in an anti- 
clockwise direction. 

When evaluating these integrals the following formula 
is useful 



1 



1 



z^'ds = 



tP- 



2nijr,sP" " T{p) 
The Lyapunov exponent is thus found to be 

ai 

For the diffusion case we assume that 

l/Ci3(z,/?) = [ai(l-z) + 0((l-z)^)] 

+/32[co + 0(l-z)] + 0(/34) 

We now need to evaluate 



(15) 



(16) 



(17) 



"^"-n + 0(n2-T) (18) 



ai 



and 



D 



Co 



(19) 



To obtain the a, b and c coefficients of this section we need 
to expand the zeta function (^ in powers of z around z ^ 
which will be discussed in the next section (HE) Then 
we will resum the series around z = 1 in section II C . 



|A| 



|A| 



j=o \j=o J 

resulting in two power series. Similarly for the diffusion 
calculation we expand 



i/CDiz) = Y[(i- 



. / ^np r"prr z"P(T 

^rip _3 \ 



6|A| 



(21) 




We restrict our attention to systems with no net drift, 
that is 



lim i(/"(xo)-io) = 



(22) 



Therefore only even powers of (3 appear in Eq. (|2l|). 

The set of coefficients we obtain in this way depends on 
the truncation used in the expansion of the infinite prod- 
uct. For truncation by topological length, we count cycles 
up to a given length TVtop . For maps with a few branches 
this number is limited to roughly of order ^ 10^, due to 
the exponential growth in the number of cycles with the 
topological length. All combinations of cycles with total 
length less than or equal to A^top are also included, as 
these contribute to the first A^top coefficients in each se- 
ries. Thus we obtain A^top exact coefficients in each series 
by topological length truncation. 
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For truncation by stability, we count cycles up to a 
given stability Amax, and combinations where the prod- 
uct of stabilities is less than A^ax- They have lengths up 
to A'ljiax, and so contribute to all of the first TVmax coef- 
ficients in each series, but are not the only contributions 
to such coefficients. They give us an approximation to 
the zeta function which for the intermittent case is more 
accurate than that obtained from the length truncation, 
but the values of the coefficients themselves are not exact 
beyond some iVoxact(Ainax), a quantity growing logarith- 
mically: iVoxact ~ logAinax- For intermittent maps, as 
the one we will consider, A^max increases as a power of 

Amax and A^max > A^exact- 

Often it is found that stability ordered cycle expan- 
sions lead to noisy results as a function of Amax- This is 
due to the breaking of shadowing pairs. For example a 
cycle AB usually gives a contribution roughly equal to 
and of the opposite sign as the combination of cycles A 
and B (we will refer to such a combination as a pseudo- 
cycle). This means the total contribution is quite small. 
The phenomenon is called shadowing, and is the main 
mechanism for the rapid convergence of cycle expansions 
in hyperbolic systems. It is still present to some degree 
in intermittent systems. However, if one such term is 
included but the other is excluded because they lie on 
opposite sides of Amax, there may be a substantial error 
generated. 

Partial shadowing which may be present can be (par- 
tially) restored by smoothing the stability ordered cy- 
cle expansions by replacing each term with inverse pseu- 
docycle stability A^^ = (Ap^---ApJ~i by S'(A)A"^ 
Here, S{A) is a monotonically decreasing function, with 
5(0) = 1 and S{A > A^ax) = 0. 

A typical "shadowing error" induced by the cutoff is 
due to two pseudocycles of stability A separated by AA, 
and whose contribution is of opposite signs. Ignoring 
possible weighting factors the magnitude of the resulting 
term is of order A"^ - (A -I- AA)"^ w AA/A^. With 
smoothing there is an extra term of the form 5'(A)AA/A, 
which we want to minimize. A reasonable guess might be 
to keep 5"(A)/A constant and as small as possible, that 
is 



5(A) 



1 - 



A 



A„ 



e(A„ 



A) 



This function still contains a non-analytic point at A = 
Amax, however the discontinuity is now in the derivative, 
not in the original function, so a smoothing error esti- 
mated by 5"(A)/A (A < Amax) is finite. We use this 
smoothing function below when evaluating the zeta coef- 
ficients, and demonstrate the improvement numerically. 



according to section 2.1, what we need are coefficients 
from some kind of (resummed) series around z = 1. We 
now describe a method of obtaining such a series, along 
the lines of Ref. fl. 



Suppose for a moment that the series X^i^o 



' has 

a radius of convergence exceeding unity. In a practical 
calculation we have only a finite number n (say, A^top 
or A^max) of coefficients Ui at our disposal. We assume 
them to be exact, the treatment of the approximate co- 
efficients from stability ordered expansions are discussed 
in sec. HI C. Then we can in principle expand it into an- 



other truncated (resummed) Taylor scries around z = 1. 



i=0 i=0 



(23) 



This leads to a linear systems of equations which is triv- 
ially invertible 



E 
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In this way one obtains the standard formulae 

y^( — X)fc log Ai + ...+log Afc 



A = 



|Ai...Afc| 



^1 |Ai...Afc| 

V^/ 1 \fc (q"l + ---+a'fc)^ 
1 Z^l ^1 IA,...Ai. 



2 E(-i)''^i^vHfr 



(24) 



(25) 



(26) 



where the sums run over all distinct pseudocycles. 

This approach is particularly cumbersome for intermit- 
tent systems where en (as well as 6, and Ci) decays ac- 
cording to some power law. Then the coefficients either 
diverges or converges slowly as n — s- oo. So, for inter- 
mittent systems the resummed series cannot be a Taylor 
series, it has to be some generalized power series. 

Assume that the asymptotic behavior of the coeffi- 
cients is a power law 



ai ^ n 



-(7-1-1) 



(27) 



Then the leading singularity is of the form (1 — z^ , and 
the simplest possible expansion would be 



^a,(l-z)^ + (l-z)'^^a,(l-z)^ 

oo 

^a^z^ (28) 



1=0 



Having only a finite number n of coefficients a; we pro- 
pose the following resummation B 



C. Resumming zeta functions 

The result of the cycle expansions in sec 2.2 is a set 
of power series of the form J2i ^i^' around z — 0. And, 



£a,(l-z)' + (l-z)'^£a,(l 

n 

^a,z^ + 0(z"+i) 



i=0 



(29) 



i=0 
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li Ua + fia + 2 = n + 1 we just get a linear system of 
equations to solve in order to to determine the coefficients 
Ui and fia from the coefficients ai. It also natural to 
require that \na + j — na\ < 1- 

The basic philosophy is to build in as much as informa- 
tion as possible into the ansatz. If the original power se- 
ries correspond to the unweighted zeta function we know 
that ao = 0. The ansatz is thus accordingly modified, we 
fix ao = and modify Ua or fia so we still get a solvable 
system of equations. 



/(i) 




■1.5 -1 -0.5 0.5 1 1.5 2 



FIG. 1. The map for a = 0.7 



III. NUMERICAL STUDIES OF AN 
INTERMITTENT DIFFUSIVE MAP 

A. The map 

In the interval x G [—1/2, 1/2), which we call the el- 
ementary cell, our model map, following Ref. jlj] takes 
the form 



f{x)^x{l + 2\2xn 



(30) 



The parameter range we consider here is a G (0, 1), where 
the Lyapunov exponent and diffusion coefficient are both 
nonzero. For any value of a, this maps the interval x G 
[-1/2,1/2) monotonically to [-3/2,3/2). Outside the 
elementary cell, the map is defined to have a discrete 
translational symmetry, 

f{x + n) — f{x) + n n G Z . 

See Fig. |l|. A typical initial x in the elementary cell 
diffuses, wandering over the real line. The map is par- 
ity symmetric, /(— i) = — /(i), so the average value of 



x„_i_i —Xn is zero, and there is no mean drift, as expressed 
in (II). 

We now restrict the dynamics to the elementary cell, 
that is, we define 

X = X - [x + 1/2] , 

where [z] is the greatest integer less than or equal to z, so 
that X is restricted to the range [—1/2, 1/2). The reduced 
map is 



fix) = fix) - [fix) + 1/2] 



(31) 



As discussed in Ref. , the intermittency of this map 
appears in the form of long cycles near the marginal point 
with power law stabilities. This is in contrast to Axiom- 
A systems for which A may be bounded by exponentials 
of the topological length. 

The map has three complete branches in the elemen- 
tary cell. Symbolic dynamics is introduced by labeling 
the branches { — ,0,-1-}. Due to the completeness of the 
symbolic dynamics the zeta functions are approximated 
by M 



and 

yCniz) 



^ ~ X! lA |i-/3 ~ ^ lA |i- 



-,n+l 



|A, 



n+l 



(32) 



(33) 



This approximation may seem crude. For instance, the 
zeta functions (32, 33) fail to preserve flow conservation. 
However in Q| we presented evidence that they capture 
the leading singularity structure correctly. This was ob- 
tained by comparing coefficients of the piecewise linear 
approximation of the intermittent map (sharing the sin- 
gularity structure with the approximation above) by the 
exact cycle expansion. The asymptotic behavior of the 
fundamental cycles is given by 



A- 



0" 



A- 



-0" 



l+l/a 



(34) 




see eg. || 4|,|14|| f or a derivation. We obtain immediately 
from d3S 



a,, 

Cn 



-1-1/" logn 



(35) 
(36) 
(37) 



This leads to the forms ( |l3|Jl7| ) with 7 = l/a as long as 
a < 1. 

For a general orbit we can only bound the stability in 
the range 

Cn^+i/" < |Ap| < (maxl/'l)"'' = (3 + 2a)"'' (38) 
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so when using stability cutoff we get for t he parameters 



and 



and Nr- 



discussed in section II E : 



A^max 



iVn- 



> 



log An 



log(3 + 2a) 



(39) 



(40) 



B. Resumming topologically ordered cycle 
expansions 

We will most of the time concentrate on the diffusion 
coefficient, but a similar analysis holds for the Lyapunov 
exponent, to which we return at the very end. 

We calculated the diffusion coefficient from resummed 
cycle expansions obtained using topological ordering as 
described in Sect. [IC with the number of coefhcients n 
determined by the maximum topological length, up to 
10. We also used the direct formula (P6[), and performed 
direct simulations with roughly the same amount of com- 
puter time. The results are shown in Fig. ^, showing that 
the resummation gives much improvement, and is consis- 
tent with direct simulation. 




FIG. 2. The diffusion coefficient at a = 0.7, from direct 
simulation (solid line), topological ordered cycle expansions 
with (diamonds) and without (plusses) resummation. 




FIG. 3. The magnitude of the zeta coefficients for stability 
ordered cycle expansions at a = 0.7 unsmoothed (points) and 
smoothed (solid line). 



The next issue is how to best make use of the infor- 
mation contained in the a„ coefficients. As pointed out 



in Sect. P B| these coefficients are not exact, but they 
give a better representation of the zeta function than the 
limited number of exact coefficients obtained from topo- 
logical ordering. In order to match the series at z = 1, 
we must again solve a linear set of equations, but the 
number of coefficients (A'max) for intermittent systems is 
much larger than for the topological ordering. We can- 
not match such a large number of coefficients in both se- 
ries, because the solution would be unstable to the errors 
in the coefficients, so we must represent the information 
contained in the a„ in the (fewer) number of degrees of 
freedom that the expansion really contains. 

There may be more than one solution to this problem; 
the solution we use here is to perform two resummations, 
the first from z = to an intermediate < z < 1, and 
the second from z = z to z = 1. 



A'max 

E 

i=0 



N„ 



which can be explicitly inverted 



iVmax 

E 



(41) 



(42) 



C. Resumming stability ordered cycle expansions 

Now we come to the central part of our numerical 
work: the resummation of stability ordered cycle expan- 
sions. First we calculate the as described in Sec. II B . 
The coefficients are all negative except So = 1, and 
their magnitudes are plotted in Fig. where we have 
used A,„ax = 10^ which corresponds to A'^cxact — 8 and 



act for n < iVcxact- The smoothed are not exact but are 
still quite accurate. For n > A'oxact we clearly see how 
the unsmoothed begins to oscillate in an irregular fashion 
where as the smoothed ones are stable for much larger n. 



With z suitably chosen, we have thus used the informa- 
tion available in the a„ approximately in proportion to 
their reliability. That is, the accurate low order coeffi- 
cients appear with large weights in the first few a„, while 
the less accurate high order coefficients appear with small 
weights. As we will see, this approach is better than one 
which simply ignores the higher order coefficients (this 
corresponds to putting z = below). 

As for the topological length truncation, the resumma- 
tion from z = z to z = 1 leads to a set of linear equations 
obtained by equating coefficients in 



E«'(i 

i=l 



zY + 



(l-z)^£a.(l 



i=0 
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(43) 



Again, we adjust and fia so as to obtain a consistent 
series in powers of z — 1 and a consistent set of linear 
equations. 



D 



D 



D 0.127 




FIG. 4. Diffusion coefficient calculated using resummed 
stability ordered cycle expansions, showing dependence on n 
(labels on curves) and z (horizontal axis). Note the scales on 
the vertical axes. 
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FIG. 5. Diffusion coefficient calculated using resummed 
stability ordered cycle expansions, showing dependence on the 
stability cutoff. The horizontal lines are the range indicated 
by direct simulation (see also discussion in text). 

We have two parameters in the double resummation 
scheme, n and z . The idea is to calculate A or D for a 
range of both parameters and look for the most consis- 
tent solution. The general behavior is shown in Fig. H. 



For each stabiUty cutoff A„iax, small values of n lead 
to a variation of D with z which has a single maximum. 
Larger values of n lead to functions that are either mono- 
tonically decreasing or oscillatory. We estimate the dif- 
fusion coefficient by finding the maximum for the largest 
value of n before monotonic or oscillatory behavior sets 
in. The convergence of this method with the stability 
cutoff is shown in Fig. 

This figure also contains the direct simulation results, 
obtained by estimating the left side of Eq. (H) for 3 x 10^ 
iterations over a sample of 3 x 10'^ trajectories, similar 
to the computer time required to find the cycles with 
A < 10'"'. The errors were obtained by looking at the 
scatter in this statistical sample of trajectories; for in- 
termittent maps the diffusion coefficient always tends to 
be too high because long intermittent episodes are not 
sampled sufRciently. Close to the phase transition at 
a = I convergence is practically logarithmic in the num- 
ber of iterations, with exponentially long times required 
to achieve convergence. For example, with the numerical 
procedure described above we find D = 0.0524 ± 0.0005 
at a = 1 where we know D = 0. Even at a = 0.7, a rea- 
sonable distance from the transition, the resummed cycle 
expansion result is more accurate than direct simulation. 

Our final value for the diffusion coefhcient at a = 0.7 
is D = 0.1267 ± 0.0003 with the resummation method. 
It is then quite compa tible with the topological order- 
discussed 



mg 
D 



m Sect 
1.262 ± 0.0003. 



[II B 



which yielded the result 
In this example, it is clear that 
resummed cycle expansions, whether ordered by topo- 
logical length or stability provide an accurate method 
of analyzing intermittent systems. Stability ordering is 
most important in more complicated systems where topo- 
logical ordering is not a realistic alternative. There is an 
additional advantage with the stability ordered expan- 
sion in that it provides a large number of approximate 
coefficients, thus facilitating a numerical estimate of the 
power law if it is not known analytically. Recall that this 
power law is used in the resummation ansatz, and was 



absolutely essential for the good result in Sect III B 



D. The phase transition 

Having gained confidence in the resummation method 
for a = 0.7, far from the phase transition at a = 1, we 
now vary a, including values for which direct simulation 
is totally impractical, due to logarithmically slow conver- 
gence. At a = 0.99 we obtain Fig. ^ for the diffusion coef- 
ficient, showing a consistent value of 13 = 0.0066±0.0001. 
Plotting D vs a (Fig. 0) we find a linear dependence near 
the phase transition at a = 1. 

Finally we performed the same analysis for the Lya- 
punov exponent, which has a similar dependence on a, 
shown in Fig. || 
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D 




series, ^ ■ Oi exp(— sZj), where the lengths of the pseudo 
orbits li are not restricted to integer values. With an 
additional resummation step at s' (corresponding to z' 
in this paper), the zeta function may be represented as a 
standard power series, thus allowing it to be matched to 
a generalized power series at s = 0. For intermittent sys- 
tems s = is again a branchpoint, and information about 
it can be obtained from the methods described in [p^ 17 
or numerically from the stability ordered expansion. 



FIG. 6. Diffusion coefficient calculated using resummed 
stability ordered cycle expansions for a = 0.99. 
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FIG. 7. The diffusion coefficient as a function of the pa- 
rameter Q, showing the approach to the phase transition at 
a = 1, beyond which D = 0. 




FIG. 8. The Lyapunov exponent as a function of the pa- 
rameter Q, showing the approach to the phase transition at 
a = 1, beyond which D = 0. 



IV. CONCLUSION 

We have demonstrated that resummed stability or- 
dered cycle expansions can provide accurate estimates 
of dynamical averages for intermittent maps, even close 
to a phase transition. This analysis could equally apply 
to maps with uncontrolled symbolic dynamics, as long as 
a reliable method exists for locating the cycles. 

Our methods can also be applied without much mod- 
ification to flows. Then the variable z is replaced by 
exp(— s) and the cycle expansion is actually a Dirichlet 
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